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The high-temperature properties of the Al(lll) surface are studied by molecular-dynamics simula- 
tion. This surface does not melt below the bulk melting point, but can be superheated. Superheating 
of metal surfaces has been recently observed in several experiments. A molecular-dynamics study 
of the structural properties reveals how after going through the superheating regime melting occurs 
over the whole crystal in a narrow temperature range. The temperature dependence of the sur- 
^ face stress, the mean-square vibrational amplitudes and the anomalous outward expansion of the 

. distance between two top layers are calculated. A transition from superheated to liquid state is 

' analyzed using kinetic description for the formation of liquid nuclei by the Fokker-Planck equation 

and conservation of heat at the liquid-solid interface. 

'Jj ■ 68.35.Rh,64.60.My,64.70.Dv 

I. INTRODUCTION 

d The question of how melting occurs, in spite of its long history and importance, remains to a large extent unanswered. 

^ Although various metastable phases of different materials are present in nature, for some time it was a common belief 
I \ that metastability is not possible at the solid-to-liquid transition. The argument was that a free surface of any material 
' O in contact with vapor acts as a nucleus of a liquid phase and suppresses superheating, i.e., nonmelting above the bulk 
^ _ melting temperature. Nevertheless, superheating was recently observed in several materials. For metal crystals it is 
especially difficult to achieve superheating. The viscosity of liquid metals is low and the liquid-solid interface rapidly 
propagates into the bulk. Several ways to suppress melting of metals were found. One way was to enclose the metal 
crystal into another metal with higher melting temperature. For example, superheating up to 62 K was observed for 
small Pb precipitates in AllJl. It was checked that the Eb particles of mean size ~ 200 A remain superheated at 18 
K above the bulk melting point for mare than 21 hoursl. The Ag spheres were coated with Au and superheated by 
25 K for a period of about one minutea. Such a behavior was modeled in molecular-dynamics (MD) simulation using 
\ the Lennard- Jones potentials and different strength of interactiond. 

. Techniques to superheat clean metal surfaces were also developed. It is known that at the temperatures below the 
' bulk melting point some metal surfaces exhibit surface melting. Using medium-energy ion scattering it was found 
, that the Pb(llO) surface melts and that the thickness of liquid layer diverges when the temperature approaches the 
' bulk melting pointo. On the contrary, Pb(lll) docs not melt and PbQOO) exhibits incomplete melting characterized 
by the presence of liquid diffusion only in one or two topmost layersu. Similar behavior show low-index surfaces of 
other fee metals, in particular surfaces of aluminunJj. It was also found that some close-packed metal surfaces under 
special conditions *nay exhibit superheating. Herman and coworkers have found superheting by 120 K for Pb(lll)B, 
90 K for Bi(0001)Q, and even 15 K for incompletely melted Pb(100)B. Superheating of these surfaces was studied 
(~| ' by time-resolved reflection high-energy electron diffraction. Melting was prevented by rapid heating with a pulsed 
O ■ laser beam. Using laser pulses of the width ~ 10^ ps it was possible to bypass melting by inducing large heating 
O ■ and co^lgg rates of about 10^^ K/s. By applying the same procedure it was not possible to ^uperheat the Pb(llO) 
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surfac Superheating was also found in MD simulation of laser-pulse irradiation for Cu(lll)0. For a development 
of new methods in materials processing by the laser beams it is important to understand the properties of superheated 
surfaces. 

I Small crystallites bounded by nonmelting facets show superheating. For example, nonequilibrium octahedral lead 
" " ' crystallites on_eraphite (made up of the (111) facets and small round parts) exhibit superheating by seveml K and for 
several hauralj. Similaii-superheating was also observed in MD high-temperature studies for Au(lll)E^, A1(111)EJ, 
Cu(lll)tj, and Pb(lll)E3. DiToUa and coworkers have performed MD simulation for deposition of a liquid aluminum 
cluster on the melting Al(llO) and nonmelting Al(lll) surfacealil. They used these results and themiodynamic 
arguments to connect superheating with the wetting angle and the non-melting induced faceting angletj. The type 
of superheating related to fcc(lll) is enabled by the exclusive presence of these non-melted surfaces on the specially 
prepared crystallites and on MD slabs. On equilibrium crystals the surface-melted facets (such as fcc(llO)) are also 
present. They act as a nucleus of the liquid phase and prevent superheating of the whole crystal. 
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Till now the studies of superheating phenomena have been mainly concentrated on the means to achieve the 
superheated state. In this work MD simulation method was used to study the properties of this metastable state, 
in particular the superheated Al(lll) surface. The analysis of the mean-square displacements, surface relaxation, 
MD particle trajectories, and surface stress has been carried out. The results of kineticj-theory based on the Fokker- 
Planck equationEZi and an analysis of the heat transport at the liquid-solid interfacet§Ej were used to estimate the 
maximum superheating temperature for aluminum. This analysis employed several quantities that were experimentally 
determined and deduced from MD simulations. In the following the MD simulation method is described in Sec. II. 
Results of simulation and discussion are presented in Subsec. Ill A. Subsection III B deals with the kinetic analysis 
of transition from the superheated to liquid state. A summary and conclusions are given in Sec. IV. 



II. MOLECULAR-DYNAMICS COMPUTATIONS 



The structure of the aluminum surface at high temperatures-.was studied by MD simulation. The interatomic 
interactions were derived fcom a classical many-body potentialE3. This form of potentials gives a proper physical 
picture of metallic bondin^j. The optimal set of parameters in the potential was founcLtty the force-matching fitting 
to ab initio electronic structure calculations using a numerical optimization procedurecEl. The data used for fitting 
were generated from different geometries, both at T = K and for finite temperatures. As a result the potential is 
characterized by a good transferability. This potential was alrcachi,|Uaed in MD simulations and it reproduced well 
experimental results for bulk and surface properties of aluminumtjE3E3. The melting point — 939 ± 5 K is in 
good agreement with the experimental value of 933.52 Ko. The calculated bulk melting temperature was precisely 
determined by simulating coexisting liquid and solid phases under constant energynJ'Ej. 

The simulation of the high-temperature properties of the Al(lll) surface started at T = K. The MD box of 1600 
particles was used. These atoms were arranged in the usual single slab geometry with the thickness of = 16(111) 
layers (i.e., ~ 35 A), and a 10 x 10 square lattice in each layer. As in other classical MD studies of surfaces, the slab 
was extended along x and y by using periodic boundary conditions and no boundary conditions were used along z. 
The three bottom layers of the slab were kept fixed to simulate the bulk. The lattice constant was changed with 
temperature according to the expansion coefficient found in MD simulation under zero pressure. The time step of 
2.64 X 10""'^^ s was used. Most runs were carried out at the constant temperature (i.e., in the canonical ensemble). 
The temperature was controlled by rescaling the particle velocities at each time step. At each temperature at least 
lO** time steps were performed to ensure thermal equilibration. In the superheating regime the length of the runs was 
much longer, i.e., up to 10^ time steps. In this regime simulations at fixed energy were also performed and the same 
results as in the canonical ensemble were obtained. 

The stress analysis often shows important changes of surface properties. The surface stress tensor atj is given by 
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where 7 is the surface free energy jper unit area, Sij is the Kronecker symbol, Cij is the surface strain tensor, and i, j 
are directions in the surface planccS. The surface stress can be obtained from the components of the pressure tensor 
calculated during MD simulationE3. In the computation of the stress MD boxes equilibrated for at least 10'^ time 
steps were used. 



III. RESULTS AND DISCUSSION 



A. Structure 



The studies of the various structural properties of the Al(lll) surface (such as the density, the static structure 
factor, the orientational order parameter, and MD particle trajectories) were done. The detailed analysis of these 
properties between K and 1500 K shows that at ^ 1120 K melting starts on a typical MD simulation time scale. 
The temperature of 1050 K was selected for a detailed analysis of the superheating regime. At this temperature the 
surface stays superheated after 10^ time steps. The static structure factor, the density and order parameter plots 
(not presented here), as well as the MD particle trajectories in Fig. |l| show that the superheated Al(lll) surface at 
1050 K and after 10^ time steps is crystalline and well ordered. It was found that in the superheating regime single 
adatoms often appear above the surface (see Fig. |l| (b)). These adatoms have a very short life time (~ 10"^^ s). 
This is in contrast with the results obtained using the same type of potentials for the high-temperature behavior of 



2 



fcc(llO) (surface melting)0 and fcc(lOO) (incomplete melting)ii. For these surfaces the formation of adatoms was not 
observed. The Al(lll) surface at T=1120 K stays superheated for a long time of 4 x 10^ MD steps. In the interval of 
(40 — 45) X 10'' time steps melting begins and develops. It was found that at this temperature melting proceeds with 
different speed in different MD runs. For example, while in one run eleven layers were melted after 42 x 10^ MD steps, 
in another one after 45 x lO** time steps only five top layers of MD box were melted. Internal energy as a function of 
temperature is shown in Fig. |^. The jump region on the caloric curve between the crystalline and liquid states is less 
abrupt than usually for melting transitions. The maximum superheating temperature estimated in this simulation is 
180 K_DiTolla and coworkers reported ^ 150 K for the same surface, using the same potential and similar simulation 
timesllJ. This is a result of different shape and size of MD boxes used in these simulations. A dependence of the 
calculated. bulk melting temperature on the size, shape, and orientation of the samples is well known problem in MD 
simulatior£J'Ej. 

The surface stress was calculated and its temperature dependence is presented in Fig. ||. The (111) surface is 
isotropic, the xx and yy components of the stress tensor are approximately equal and therefore the stress is represented 
by their mean values. The stress of 0.056 eV was found for the relaxed surface at T = K. TIik is 3.3 times less 
than the average stress for Au(llO) at the same temperature and for the same type of potentialEd. The calculated 
stress for Al(lll) at T = K is in excellept agreement with 0.059 eV A"'^ obtained for the same surface in MD 
simulation using the Sutton-Chen potentiaS. These MD results should be compared with twp-|a& initio electronic 
structure calculations for Al(lll). The value of 0.078 eV was found by Needs and GodfreytHi and 0.090 eV A~^ 
by FeibelmanEj. Schmid and coworkers calculatcjd, surface stress for some low- index fee metal surfaces using MD 
simulation and effective medium theory potentialsEa. They found that MD generally gives lower values compared to 
electronic structure calculations. Figure 3 shows that surface stress for Al(lll) is almost constant between T = K 
and T = 900 K. The stress is always tensile. The components of any type of stress (or pressure) tensor calculated in 
MD simulation exhibit large fluctuationsES. In Fig. 3 the errors increase with the temperature and the maximum of 
statistical uncertainty for the points is 20%. Same trend in the distribution of the values (i.e., that the stress remains 
almost constant over the wide temperature region and also that pronounced scatter of the data exists) was found in 
MD simulation for the (111) surface of the Lennard- Jones crystalEII. It is well known that surface stress [see Eq. (1)] 
for a liquid becomes equal to the surface free energy (i.e., d^/dtij — 0). The value 0.05 eV A~^, obtained in this 
simulation for the liquid surface at T = 1500 K, is equal to the surface free energy of liquid aluminuirO. 

The mean-square vibrational amplitudes for the surface atoms are shown in Table |. The values of mean-square 
displacements are presented at two temperatures below the bulk melting point (300 K and 900 K) and also at 1050 K, 
i.e., for a typical temperature in the superheating region. The last row of this table is obtained for the solid surface 
at the onset of melting at T = 1120 K. The valips-fottnd for temperatures below the bulk melting point are in a 
good agreement with other MD simulation resultsEjCJ^Lll. The mean-square vibrational amplitudes at the end of the 
superheating region (not studied elsewhere) are ^ 10 times larger than below the bulk melting temperature. The 
mean-square amplitudes of vibration are isotropic at all temperatures. Most MD simulations of the highptemperature 
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surfaces also show iiiat the mean-square vibrational amplitudes are isotropic! In recent 
and MD simulationEfl of two fcc(lOO) surfaces it was found that the mean-square amplitudes 
of vibration are anisotropic. Out-of-plane vibrational amplitudes were found to be smaller than in-plane ones. This 
reveals lateral disordering typical for incomplete melting. Isotropic mean-square vibrational amplitudes found here 
for Al(lll) show that after superheating regime melting proceeds in similar way along all directions. 

Most low-index metal surfaces relax inwards at lower temperatures. It is known that the Al(lll) surface exhibits 
an anomaly in relaxation, i.e., tkat at low temperatures the distance between the two outmost layers expands in 
comparison with the bulk layer£j. The potential for aluminum used in this paper gives good . desc ription of the 
Al(lll) relaxation, i.e., it gives the experimental value for surface relaxation -1-0.9% at T = 90 KE3Ej. In this work 
the temperature dependence of surface relaxation was studied and results are also presented in Table ^ It was found 
that this unusual outward expansion of the distance between the two top layers increases with temperature. This 
increase is ~ 1% along the superheating temperature region. The maximum of surface relaxation is -|-3.3% at the 
onset of melting. 



B. Metastability 

The fundamental questions of conditions and limits for the existence of any metastable state deserve further in- 
vestigation. For metastable Al(lll) it is important to consider a maximum of the superheating temperature as a 
function of aluminum properties. For fee metal crystals the maximum of superheating should occur at the close packed 
non- melted (111) surfaces. MD simulation presented above shows that the mean-square vibrational amplitudes are 
isotropic at the transition from the superheated to liquid state (see Table |). Moreover, drops form everywhere, and 



3 



the liquid front rapidly propagates into the bulk. In the following estimate of the maximum superheating temperature 
small anisotropy of the solid-liquid interface free energy is not taken into account. Therefore, the spherical liquid 
nuclei are analyzed. 

The kinetie, theory of first order phase transitions provides a description of transformation from a metastable to 
stable phasdm. This transition proceeds via fluctuation induced formation of the nuclei of the stable phase. If the 
radius r of the nucleus is smaller than some critical value Vc such nucleus disappears, if it is bigger the nucleus grows. 
The radius of the spherical liquid critical nucleus is 

2a 

L T-Tm ^ ' 

where a is the solid-liquid interface free energy and L is the latent heat of meltin^^. It is possible to describe the 
growth of the nucleus by the kinetic Fokker-Planck equation for the distribution function f{r,t) 

df _ _dw 

dt~ 9r ' 

where w is the fiux density. The phase transition corresponds to the stationary solution of Eq. (3), where w = const. 
Such solution isEJ 

/ rv 

w^2^-B{r,)h{T,). (4) 

In this equation B is diffusion type coefficient and fo{rc) is 

/o(r-c) = const exp ^^^^^ . (5) 

Equation (4) gives the rate for formation of a critical nucleus. The goal here is to calculate the maximum of super- 
heating, i.e., the temperature at which melting occurs. Therefore, at this temperature w ~ 1. The coefficient B(rc) 
in Eq. (4), is given byllZl 

„, ^ T dr 

87ra(r — rc) at 

and can be calculated [for r — > Tc as in Eq. (4)] by considering a particular metastable state. 

For a transition from ±h£_|Superheated to liquid state it is necessary to consider the diffusive transport of heat at 
the liquid-solid interfacetSEj. It is convenient to define the dimensionless field 

« = ^(Too-r), (7) 

where Too is the temperature of the solid infinitely far from the growing nucleus and Cp is the specific heat. Then the 
diffusion equation is 

where D is the thermal diffusion constant at T ~ r,„. One boundary condition is the heat conservation at the 
liquid-solid interface 

Vn = -Dfi ■ vw, (9) 

where n is the unit normal directed outward from the nucleus and Vn is the normal growth velocity. There is also a 
requirement of local thermodynamic equilibrium that gives the dimensionless temperature ui at the interface 

Ui = A - dnK,, (10) 

where 

A = ^{T^-Tra) (11) 
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is proportional to superheating. In the second term in Eq. (10) (i.e., in the Gibbs-Thomson correction for the melting 
temperature at a cwved surface) K. is the curvature and c?o = oiCpTmL^^ . The solution for the growth rate of a 
spherical nucleus ist3 

-^(Ar,-2do). (12) 

Therefore, using Eqs. (5), (6), (12), and Too — T, condition w ~ 1 for Eq. (4) gives 

aT^/\T -T,r,f ei.Y'i-lh/TiT -T^fW = 1, (13) 

where a = DLCp/lQira^/^T^ ajid 6 = leTra^T^/SL^. Equation (13) was solved numerically using the experimental 
data for interface free, energy da, latent heat of melting Zo, bulk melting temperature TmEj, specific heat CpE3, and 
diffusion constant £cd. The value ~ 23 K was obtained for the maximum of superheating. This value is smaller than 
the maximum superheating temperature of 180 K estimated in MD simulation for a model of the Al(lll) surface. 
Part of this disagreement is the result of the approximate kinetic analysis and even possibly poor accuracy of some 
experimental data for Al. There is also a possibility that the value for the maximum of superheating found in MD 
simulations is in part the result of a limited time evolution. Although long runs of 10^ time steps were performed in the 
superheating regime, much longer (nowadays not feasible) simulation times may give the maximum of the superheating 
temperature in accordance with the kinetic result. Using tl:te-|kinetic analysis presented above, the same literature 
sources for a, L, Tm, Cp, and diffusion constant D from Rcf.ca, the value of ~ 47 K was obtained for the maximum 
superheating temperature of Pb_For aluminum a similar Fokker-Planck analysis was done using quantities deduced 
from the simulations. As in Ref.Ell, it was found that the bulk energy is slightly nonlinear function of temperature and 
that therefore it can be fitted by the second-order polynomial. Using this procedureEl, the value of Cp = 1120 J/(K 
kg) was obtained, whereas the experimental value for the specific heat is 902 J/(K kg) [23]. For diffusion constant 
the value of D = 0.3 x 10~^ err? j s was found (see alsoL5). This is much lower than the experimental vakie 3 x 10~^ 
cmLjxA. When in Eq. (ISLthe quantities obtained in MD simulations for Cp, D, L — 0.105 eV/atomE3, = 939 
KEHI, and a = 10 meyA^^ocj were used, then 48 K was calculated for the maximum superheating temperature of 
aluminum. It is important to point ai:d_Lhat smaller values for the maximum superheating temperatuce-were found 
in the experiments on the crystalliteaj'Erllj, whereas larger superheating was induced by a laser beanJjQllj. 

IV. SUMMARY AND CONCLUSIONS 

The properties of a surface in the superheated state were studied using MD simulation and a reliable many-body 
interatomic potentiaL.|_Superheated Al(lll) was used as a model of superheated surfaces, such as ones obtained by a 
pulsed laser bearnETErE^, or on the crystallitestj. A detailed analysis of the Al(lll) surface from room temperature 
to 1500 K (i.e., well above the bulk melting point) was carried out. The results for lower temperatures are in a good 
agreement with available experimental observations and other MD simulations and electronic structure calculations 
for metal surfaces. It is possible to superheat the sample by 180 K for typical longest simulation times used in 
classical MD (> 2.5 ns). In the superheating regime the Al(lll) surface is remarkably well ordered, although single 
adatoms sometimes appear. The sample melts over the narrow temperature interval. Anomalous outward expansion 
between two top layers increases slowly with the temperature: from -1-0.9% at T = 0, up to -1-3.3% at the end of the 
superheating region. The mean-square vibrational amplitudes are isotropic for all temperatures and ~ 10 times larger 
at the end of the superheating region than below the bulk melting point. It was shown that kinetic theory based 
on the Fokker-Planck equation and analysis of heat conservation at the liquid-solid interface for aluminum gives the 
maximum superheating temperature of 23 K when experimentally determined parameters where used. The maximum 
superheating temperature of 48 K was obtained when parameters deduced from MD simulations where applied in the 
Fokker-Planck analysis. This analysis of kinetics and process of disordering observed in MD simulation shows that 
superheated Al(lll) and bulk Al below the surface are an example of metastability at the salid-to-liquid transition. 
The kinetics of this transition is the same as in other better known examples of metastabilityt3. Superheated surfaces 
of other metals should exhibit similar behavior. 
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TABLE I. Mean-square vibrational amplitudes for the Al(lll) surface (in units of A'^): (along the [110] direction), Uy 
(along the [112] direction) and ul (along the vertical axis). The surface relaxation di2 between the two top layers is also shown 
(in %). 



Temperature 


ul 


ul 


ul 


di2 


300 K 


0.014 


0.022 


0.018 


+1.1 


900 K 


0.074 


0.064 


0.075 


+2.2 


1050 K 


0.118 


0.119 


0.101 


+2.4 


1120 K 


0.839 


0.708 


0.731 


+3.3 



FIG. 1. Particle trajectories showing the superheated surface after 10® MD steps of time evolution at 1050 K (i.e., 110 K 
above the bulk melting temperature): (a) top view, (b) side view with an adatom above the surface. Trajectory plots refer to 
a time span of ~ 3 ps and include only moving atoms. 



FIG. 2. The caloric curve. The vertical line represents the calculated bulk melting temperature. 



FIG. 3. Surface stress as a function of temperature. The dashed vertical lines enclose the superheating region. 
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This figure "Figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9704045vl 
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